Number of Transition Frequencies of a System Containing an Arbitrary Number of 

Gas Bubbles 
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"Transition frequencies" of a system containing an arbitrary number of bubbles levitated in a 
liquid are discussed. Using a linear coupled-oscillator model, it is shown theoretically that when the 
system contains N bubbles of different sizes, each bubble has 2N — 1 (or less) transition frequencies 
which make the phase difference between an external sound and a bubble's pulsation n/2. Fur- 
thermore, we discuss a discrepancy appearing between the present result regarding the transition 
frequencies and existing ones for the resonance frequencies in a two-bubble case, and show that 
the transition frequency, defined as above, and the resonance frequency have a different physical 
meaning when N > 2, while they are consistent for N = 1. 



Multibubble dynamics in an acoustic field have been 
studied in a variety of fields p|-|T3|], and the mutual ra- 
diative interaction of bubbles is known to change their 
acoustic properties [Q-^|j8[||,|T^] . Shima jjj derived a the- 
oretical formula for the resonance frequencies of two mu- 
tually interacting gas bubbles levitated in an inviscid liq- 
uid, and showed that a bubble interacting with a neigh- 
boring bubble has two resonance frequencies, which re- 
veal an upward or downward shift as the bubbles ap- 
proach each other. Zabolotskaya (|] showed theoretically 
that variations in the resonance frequencies may cause 
changes of phase difference between bubbles, and may, 
in turn, sometimes give rise to reversal of the sign of the 
secondary Bjerknes force, which is an interaction force 
acting between pulsating bubbles jlO| . In the studies de- 
scribed in refs. 3 and 4, Doinikov and Zavtrak produced 
results that were almost equivalent to those of Zabolot- 
skaya by employing a mathematical model in which the 
multiple scattering between bubbles is described using 
a Legcndrc expansion function. Their results are con- 
sidered to explain the "bubble grape", which is a sta- 
ble structure formed by bubbles in an acoustic field |Q. 
Feuillade §], using a self-consistent model, showed that 
when two bubbles are pulsating in antiphase and radi- 
ation damping is dominant rather than thermal or vis- 
cous damping, the bubbles may show a superresonant re- 
sponse P| because of the reduction in the effective damp- 
ing. Ye and Alvarez J7], using a self-consistent model, 
carried out numerical experiments on acoustic wave prop- 
agation in a liquid containing randomly distributed iden- 
tical bubbles, and showed that the localization of acoustic 
waves can take place in a range of frequencies. It was con- 
sidered that the multiple scattering of sound due to the 
bubbles has a strong influence on this phenomenon. Met- 
tin et aJL Q and Doinikov || examined some influences 
of the interaction of two bubbles in a strong acoustic field 
on the radiation forces (the secondary and primary 
fbi] , fl5f Bjerknes forces) acting on the bubbles; these stud- 
ies are related to multibubble sonoluminesccncc l|lq,|l( 



The aim of this paper is to provide preliminary but sig- 
nificant discussions regarding the "transition frequency" 
of a system containing an arbitrary number of bubbles 
having different radii, by employing a simple theoretical 



model. In the present study, the transition frequencies of 
a bubble in the system are defined as the frequencies of 
an external sound for which the phase difference between 
the sound and the bubble's pulsation becomes n/2 p2[ . 
This definition is based on the well-known phenomenon 
of a single bubble case in which the phase difference be- 
tween an external sound and the bubble pulsation be- 
comes 7t/2 at its resonance frequency (more correctly, at 
its natural frequency) [ ^3|p^ ]. This study has two main 
purposes: (1) To determine how many transition frequen- 
cies a bubble in the system has. As has been reported, 
it can be considered that the phase shifts of the bub- 
ble pulsations due to the radiation coupling give rise to 
a large variety of phenomena that cannot be predicted 
by a single-bubble model (see, e.g., refs. 1-9). However, 
there exist only a few studies on the transition frequen- 
cies defined as mentioned above. Almost all of the pre- 
vious works (e.g., refs. 1-5) concern the resonance fre- 
quencies. Very recently, Ida Jl^], using a coupled oscilla- 
tor model, showed theoretically that a bubble interacting 
with a neighboring bubble has three transition frequen- 
cies. In the present paper, by extending this approach, 
we derive the number of transition frequencies of a bubble 
interacting with an arbitrary number of bubbles. (2) To 
elucidate a difference between the transition frequency 
and the resonance frequency. As suggested already, the 
transition frequency and the resonance frequency have 
different physical meanings when the number of bubbles, 
N, is larger than 1. As briefly reviewed above, Shima |l| 
and Zabolotskaya j| predicted the existence of two reso- 
nance frequencies per bubble when N = 2. This number 
of resonance frequencies is not equal to that of the tran- 
sition frequencies given by Ida [l2] ]. In the present paper, 
we show that in general, the numbers of those frequen- 
cies arc different. This result shows that when N > 2, 
the reversal of the phase of the bubble pulsation (e.g., 
from in-phase to out-of-phase with an external sound) 
may occur not only at the resonance but also at other 
frequencies, unlike the case of N = 1. 

At first, we derive a system of linear differential equa- 
tions that describes the pulsation of coupled bubbles and 
that is almost equivalent to those used in refs. 1, 2, 5, 12 
and 19. Let us assume that N bubbles are levitated in a 
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viscous liquid and the equilibrium radius of each bubble 
(assumed to be constant ^FJ) is R^ (i — 1,2, ■■■,N). 
When the amplitude of the external sound is weak, the 
bubble pulsation can be described by a linear second- 
order differential equation p|J22,0 
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pRio 



(1) 



where it is assumed that the sphericity of the bubbles 
is maintained and the time-dependent bubble radii Ri (t) 
can be represented by Ri(t) = R40 + ei(t) (\a\ <C Rio), 
and uj i0 (= ^/[3kP + (3k — l)2a/ R i0 ]/ pR 2 Q ) is the reso- 
nance (angular) frequency, Si is the damping factor p3| , 
Pdi is the driving pressure acting on bubble i, p is the 
density of the surrounding liquid, k is the polytropic ex- 
ponent of the gas inside the bubbles, Pq is the equilibrium 
pressure, a is the surface-tension coefficient, and the over 
dots denote the derivation with respect to time. When 
more than one bubble exists, the driving pressure p<a is 
given by the sum of the pressures of the external sound, 
p ex , and the scattered sound waves due to other bubbles, 
p Bi j, as pi 
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Pdi=Pcx+ ^ Psij~Pcx+ ^2 
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where the surrounding liquid is assumed to be incom- 
pressible, and Dij is the distance between the centers of 
bubbles i and j. By substituting eq. @ into eq. (|l|), one 
obtains 
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ei + u; l0 ei + Siii = — Tr'^i' ( 3 ) 

pRio R i0 j J^ j _ ti >>', 

This kind of system of differential equations is called a 
coupled-oscillator model or a self-consistent model. Now 
we assume that p ex = Pexp(iwt) and e, = ^ exp(i wi), 
where P is a real, positive constant and Ai denotes com- 
plex amplitudes. Those assumptions reduce eq. (^) to 



JV 



R i0 [(co 2 -u? )-iu5 i }A i + TT w2A - 



p 



(4) 



This system of equations can be written in a matrix form, 

MA = B, (5) 

where M is an N x N matrix whose elements, rriij, are 
defined as 

i? i0 [(w 2 - ujf ) - iujSi], for i = j, 
= <? R 2 
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j0 Lu 2 , otherwise, 



A={A 1 ,A 2l --- 1 A N ) T andB= (P/p)(l, 1, • • • , 1) T . 

Let us analyze the matrix equation. The solution of 
eq. (||) is expressed as 



A = M~ 1 B. (7) 
As is well known, M _1 can be written in the form of 

C 



M = 



Ml' 



(8) 



where |M| and C are the determinant and the cofactor 
matrix, respectively, of M. From the definitions of the 
determinant and the cofactor of a matrix, one determines 
that 



deg(|MH|) = 2N, 
deg(C(u,))=2(iV-l). 



(9) 
(10) 



Here we should note that deg (mij(u>)) = 2. Substitution 
of eq. (H) into eq. (Q) yields 



A = 



IMI 



-CB. 



As defined previously, the transition frequencies of bub- 
ble i are derived by an equation of 



Assuming that 



Re(^ z ) = 0. 

|M| =a + ib, 
gt = Ci + id, ; , 



(11) 



(12) 
(13) 



where a, b, Cj, and di are real values and gi is an element 
of CB, yields 

Cj + i d% aci + bdi + i (adi — bci) 



and 



a + i b 



Re(A, 



a 2 + b 2 



aci + bdi 
a 2 + b 2 



(14) 



Using this equation and assuming that |M| =/= 0, i.e., no 
case exists where both a = and 6 = are true (this as- 
sumption is physically valid because this guarantees the 
absence of a singular solution) , eq. (|i"l| ) is reduced to 



ac; 



bdi = 0. 



(15) 



Now we examine how many roots this equation has. From 
eqs. ([)]) and (|h]) and definitions (jlj) and (|l3|), we know 
that deg(aQ + bdi) = 2N + 2(N - 1) = 4N - 2. Fur- 
thermore, using definition 1q), we can prove that eq. (||) 
contains terms of only even orders with respect to lu as 
follows: the real parts of m^- are of the second order and 
the imaginary ones are of the first order with respect to 
w; thus, only terms of even orders remain in the real 
part of Ai, because an even-order component given by 
the imaginary parts of rriij, e.g., iuiSi x iu52 = —6i52W 2 , 
becomes a real value and an odd-order component, e.g., 



2 



\ui8\ x ibj52 x iluS?, = — i<5i(52^3k-> 3 , becomes an imaginary 
value. As a result, eq. (pJ5|) can be represented by 



F(X) = + bd t = 0, 
and one knows that 

deg(F(X)) = 2iV-l, 



(16) 



(17) 



with X = u! 2 . These equations predict that a bubble has 
2N — 1 (or less) transition frequencies. (Note that only 
a positive root is physical.) When N = 2, for example, 
2N —1 = 3; this result is in agreement with that given 
by Ida @. 

To confirm the above result, we show here a numerical 
result for N — 3. The parameters used are Rio — 1 /^m, 
i?2o = 1-5 pm, R30 = 2.5 /xm, p = 1000 kg/m 3 , k = 1.4, 
Po = 1 atm, and a = 0.0728 N/m. The damping factor 
is set to Si = 4p/pR 2 , where p (= 1.137 x 10~ 3 kg/(m 
s)) is the viscosity of water at room temperature (i.e., 
the viscous damping is adapted, which is dominant for 
small bubbles used, e.g., in medical applications |Tl| l and 
in the experiments on sonoluminescencc p^Jl3] | ) , and Dij 
is determined by — s(Rio + Rjo)- For s = 2.0, we, as 
expected, obtain five (=2x3 — 1) transition frequencies 
of bubble 1, wi = (1.054, 0.723, 0.635, 0.411, 0.328) x w w , 
while for s = 10.0, only one, w\ = 1.002ojio, is found. 
Such a dependency of the number on the distance be- 
tween bubbles can also be found in a two-bubble case 
jl2| . Detailed discussion regarding three- or more-bubble 
cases will be given in a future paper. 

We perform here a comparative study between the 
present theory and previous ones in order to clarify a 
difference between the transition frequency and the reso- 
nance frequency. Shima Jl]] and Zabolotskaya [|| derived 
the same formula for estimating the resonance frequen- 
cies in the case of N = 2, expressed as 



{X-u\ ){X-ul Q ) 



(18) 



In their mathematical models, the sphericity of bubbles 
and incompressibility of the surrounding liquid were as- 
sumed as in the present one, but the damping factor was 
neglected. This formula predicts two resonance frequen- 
cies; apparently, this result seems to contradict the result 

§5n by Ida, which predicts three transition frequencies 
. We next discuss the origin of this discrepancy. 
Generally, when no damping factor exists (Si = 0, 
Im(mjj) = 0), the resonance frequencies are determined 
by 



I Ml = 



(19) 



in the case of this study. This condition gives rise to 
the well-known singularity in the bubble pulsation at the 
resonance frequencies (infinite amplitudes and discontin- 
uous phase reversals) pE| , [L4| . This formula is reduced to 
the well-known one for a single-bubble case, u 2 —ujf = 0, 



and to eq. ( |18| ) for a double-bubble case. On the other 
hand, in the present study, the transition frequencies of 
bubble i in the case of Si = are derived by the following 
independent equations: 



(20) 



and 



Sign I lim AAu + A) ^ Sign lim AAu + A) 



(21) 



where 



Sign(/) = 



1 for / > 0, 
-1 for / < 0. 



Equation ( pl| ) means that the phase of bubble i is re- 
versed at u> satisfying this equation. (In the limit of 
Si —> 0, the condition that the phase difference is tt/2 
converges to eq. ( f2l| ) at the resonances because of the 
singularity mentioned above. Namely, at the resonances 
those conditions are essentially equivalent.) Equations 
( pp| ) and (|l]), respectively, may be reduced to gi = 
and |M| = 0, or 

\M\ gi = Q. (22) 

Furthermore, using the definitions of eqs. ( |i"2| ) and jl^), 
eq. ( p2|) can be reduced to 

aci = 0. 

This equation is in agreement with eq. ( |i"5| ) when 6 = 
and di = (these conditions are naturally satisfied be- 
cause of Im(my) = 0). The degree of eq. jlS| ) with re- 
spect to X is N as discussed above, while that of eq. ( p^ ) 
is N + (N - 1) = 2N - 1 . When N > 1 , these two degrees 
are not equal to one another. This result reveals that (1) 
"resonance frequency" and "transition frequency" have 
different physical meanings, (2) the number of the transi- 
tion frequencies is, in general, larger than that of the res- 
onance frequency, and (3) other conditions under which 
the phase reversal of bubble pulsation takes place exist in 
addition to those for the resonance. The transition fre- 
quencies given by |M| =0 cause the resonance response 
of the bubbles, while those given by gi = do not cause 
it. (The respective transition frequencies may correspond 
to the series and the parallel resonance frequencies of an 
electromagnetic circuit, which give rise to a zero and an 
infinite impedance, respectively.) 



Next, we discuss the case of D; 



00 and Si 



0. 



i.e., in which the non-diagonal elements of M and the 
imaginary parts of its diagonal elements have quite small 
absolute values. These settings result in 



N 

a= I] Rjo(u 2 ~ Lo 2 ) + si, 
3=1 



b = e2, 



P 



N 



Ci = — TT Rjo(u 2 - oj 2 ) + e3, d t = £4, 
P , 
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where el ~ e4 are real values whose absolute values are 
quite small. They reduce eq. (|i~5| ) to 
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This equation shows that for Dij — > oo, among 2A r — 1 
transition frequencies of bubble i, only one converges to 
LOiQ, while two of the remaining ones converge to lojq 
(j ^ i). This result is consistent with that for N = 2. 

Lastly, we briefly discuss a special case. Here we as- 
sume that all bubbles are identical (i?io = R2Q — ■ ■ ■ = 
Rno, Si — 82 = ■ ■ ■ = <5/v) and, furthermore, that the 
distances between all pairs of bubbles are uniform. (The 
latter assumption may be realizable up to N — 4.) When 
N = 3, this setting corresponds to the "mode A" dis- 
cussed in Sec. Ill of ref. 5. Based on these assumptions, 
the elements of matrix IMI become 



physical meaning of the existence of an odd number of 
transition frequencies may be explained as follows: a bub- 
ble, even if interacting with neighboring bubbles, pulsates 
in-phase or out-of-phase with an external sound if the 
driving frequency is much lower or much higher, respec- 
tively, than the resonance frequencies of the bubbles; to 
interpolate those two extreme states consistently, an odd 
number of the phase reversals is needed.) This result 
is important especially for subjects in which the phase 
change of bubbles plays an important role, such as those 
studied in refs. 2-4, 8 and 9 and investigated using only 
the resonance frequencies. (In ref. 26, the author already 
applied the theory for the transition frequency to the in- 
vestigation of the reversal of the sign of the secondary 
Bjcrknes force, and found that the theory gives an alter- 
native interpretation of the phenomenon, which may be 
more accurate than the previous ones given in refs. 2-4.) 
We expect that more detailed mathematical discussions 
regarding the matrix equation will supply a still richer 
understanding of the physics of multibubble dynamics in 
an acoustic field and acoustic wave propagation in bubbly 
liquid. 



run = TO22 = • • • = Wjvjv = i?io[(^ 2 - wf ) - iuSi], 
mi2 = mn = TO13 = • ■ • rriNN-i 



K> 2 
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Namely, both the diagonal and the non-diagonal elements 
become uniform. As a result, one obtains A± = A2 = 
■ • • = An and 



Ax 



th\\ + (N — 1)toi 



Only one transition frequency (corresponding to the res- 
onance frequency) is given by this formula, although 
2N — 1 ones are predicted in general. This result points 
out that the numbers of both the resonance and the tran- 
sition frequencies are reduced under certain conditions 
and that such an excessive idealization prevents the accu- 
rate understanding of multibubble dynamics in an acous- 
tic field. 

In summary, we have studied the transition frequen- 
cies of a system containing an arbitrary number (N) of 
bubbles of different sizes by employing a linear coupled- 
oscillator model. In the present study, the frequency 
of an external sound in the case where the phase dif- 
ference between the external sound and a bubble's pul- 
sation becomes n/2 was called the transition frequency. 
The present theory predicts that a bubble in the system 
has 2N — 1 transition frequencies, while it has N reso- 
nance frequencies; namely, the well-known theory for a 
one-bubble case in which the phase reversal of a bubble's 
pulsation takes place around its resonance frequency is 
not absolutely true in a two- or more-bubble case. (The 
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